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ABSTRACT 

We investigate the constraints on primordial non-Gaussianity with varied bispectrum shapes 
that can be derived from the power spectrum of galaxies and clusters of galaxies detected in 
future wide field optical/near-infrared surveys. Having in mind the proposed ESA space mis- 
sion Euclid as a specific example, we combine the spatial distribution of spectroscopically 
selected galaxies with that of weak lensing selected clusters. We use the physically motivated 
halo model in order to represent the correlation function of arbitrary tracers of the Large Scale 
Structure in the Universe. As naively expected, we find that galaxies are much more effective 
in jointly constrain the level of primordial non-Gaussianity /nl and the amplitude of the mat- 
ter power spectrum <x 8 than clusters of galaxies, due to the much lower abundance of the 
latter that is not adequately compensated by the larger effect on the power spectrum. Never- 
theless, combination of the galaxy power spectrum with the cluster-galaxy cross spectrum can 
decrease the error on the determination of /nl by up to a factor of ~ 2. This decrement is par- 
ticularly evident for the less studied non-Gaussian bispectrum shapes, the so-called enfolded 
and the orthogonal ones. Setting constraints on these models can shed new light on various 
aspects of the physics of the early Universe, and it is hence of extreme importance. By com- 
bining the power spectra of clusters and galaxies with the cluster-galaxy cross spectrum we 
find constraints on primordial non-Gaussianity of the order A/nl ~ a few, competitive and 
possibly superior to future CMB experiments. 



1 INTRODUCTION 



It is now commonly accepted that the formation of structures in 
the Universe originated from seed density fluct uations in the dark 
matte r fluid that were laid down during inflation ( Hawking & Moss 
Il983h . Gravitational instability has the effect of amplifying den- 
sity perturbations, that at a given point enter the non-linear regime 
and collapse to form clumps, voids, and more in general the com- 
plex Large Scale Structure (LSS henceforth) that we observe to- 
day. For a given distribution of the initial conditions, the statis- 
tical properties of the LSS are determined uniquely by the sub- 
sequent expansion history of the Universe that is, ultimately, on 
its matter and energy content. In the past decade much effort has 
been directed toward understan ding the effect of dark energy on 
the formation of s t ructures (see[C unha, Huterer, & Frieman 2009; 
Grossi & Springell 120091; ISamushia & Ratral 120091 ; |Pe Boni et all 
2010; Mortonson, Hu, & Huterer 2010 for recent examples) in or- 



specifically about the shape of the primordial density fluctua- 
tions distribution. While the simplest models of inflation (single 
slow-rolling scalar field) predicts this distribution to be virtually 
indistinguishable from a Gaussian, a plethora of more complex 
models have been proposed that, in addition to solving the standard 
cosmological problems, allow for significant and possibly scale 
dependent deviations from Gaussianity. The study of non-Gaussian 
cosmologies and the effect they have on the formation and evolu- 
tion of cosmic structures is thus extremely important in order to 
rule out inflationary models, and hence to have a better handle on 
the physics of the early Universe. Moreover, studies on possible 
detectability of primordial non-Gaussianity are very timely, given 
the recent claims for a positive skewness of the primordial density 
fluctuations distribution coming from the Cosmic Microwave 
Background (CMB) a nd from the ang ular correlation function of 
radio-selected quasars JXia et alfeOld) . 



The problem of constraining deviations from primor- 
dial Gaussianity by means different from the CMB intrinsic 
anisotropics has recently attracted much efforts in the litera- 
ture, with st udies directed towards the ab u ndance o f non-linear 
After several pioneering works jMessina etal] Il99d ; structures ( Matarrese, Verde, & Jimenez 



der to gain insights on its nature, specifically whether it truly is a 
cosmological constant or it does have some kind of evolution with 
cosmic time. 



P« 

Moscardini etal.lll99ll ;IWeinberg & Coldll992t) . only recently has 
the question of the initial conditions gained renewed attention 
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tegrated Sachs-Wolfe effect dAfsh ordi & TolTej] I2008J : 
ICarbone. Verde. & Matarresell2008 ), Lya flux from low-de nsity in- 
tergal a ctic medium dViel et al.1 2009), 21 — cm fl uctuations dCooravl 
2006; Pillepic hTPorciani. & Matarresd |2007|) and reionization 
dCrociani et alj|2009h . 

In this work we focused attention on the spatial distribution 
of galaxies and galaxy clusters as tracers of the LSS. We were par- 
ticularly interested in comparing the performances of the power 
spectra of the two individual tracers in constraining primordial 
non-Gaussianity, and evaluate the improvements in forecasted con- 
straints given by the addition of the cluster-galaxy cross power 
spectrum. Throughout the paper we assumed a fiducial f uture all- 
sky optical/near-infrared s urvey on the model of Euclid dLaureiisI 
l2009l ; lBeaulieu et aljfeoid) . In order to fully exploit the potentials 
of both the imaging and the spectroscopy part of Euclid we consid- 
ered galaxies as spectroscopically selected according to their Ha 
flux, and galaxy clusters as selected to be the high signal-to-noise 
ratio (S/N) peaks in full-sky cosmic shear maps. This approach 
has the advantage of allowing treatment of the galaxy and clus- 
ter samples as independent. The results obtained here are relevant 
for other planned missions with a similar concept to Euclid, such 
as the NASA Wide Field Infrared Survey (WFIRST). 

The rest of the paper is organized as follows. In Sections [2] 
and [3] we summarize the non-Gaussian models that we have ex- 
plored in this work, as well as their effect on the mass function 
and large scale bias of dark matter halos. In Section [4] we de- 
scribe the halo model, the physically motivated framework that we 
adopted for modeling the power spectrum of clusters and galaxies 
as well as the cluster-galaxy cross spectrum. In Section [5] we sum- 
marize our results and in Section[6]we discuss them, with particu- 
lar emphasis on alternative survey configurations. Finally, in Sec- 
tion [7] we present our conclusions. Throughout this work we re- 
fer to the fiducial cosmological model as the one defined by the 
parameter set derive d by the latest analysis of the WMAP data 
dKomatsu et alj2010f) . This means that density parameters for mat- 
ter, cosmological constant, and baryons are equal to fl m = 0.272, 
£^a,o - 0.728, and fl b = 0.046, respectively, the Hubble constant 
equals h = Hq/( W0 km s _1 Mpc~') = 0.704, and the normalization 
of the matter power spectrum is set by cr 8 = 0.809. 



2 NON-GAUSSIAN COSMOLOGIES 



Exten s ions of the most standard model of inflation dStarobinskiil 
1 19791 : iGuthl Il98ll ; iLindel Il982l) can produce substantial de- 
viations from a Gaussian di stribution of primo r dial density 
and potential fluctuations (see iBartolo et alj 1 20041 : IChenl 1 20 1 Ol ; 
Desjacques & Seliak 20101 for recent reviews). The amount and 
shape of this deviation depend critically on the kind of non-standard 
inflationary model that one has in mind, as will be detailed later on. 

A particularly convenient (although not unique) way to de- 
scribe generic deviations from a Gaussian distribution consists in 
writing the gauge-invariant Bardeeen's potential $ as the sum of a 
Gaus s ian random field and a quadratic correct ion ( Salopek & Bond 
1990 ]: iGangui et al] 1 1994 IVerde et al.l l200d : iKomatsu & Spergell 
200 lL according to 



The parameter /nl in Eq. (TJ determines the amplitude of non- 
Gaussianity, and it is in general dependent on the scale. The 
symbol * denotes convolution between functions, and reduces to 
standard multiplication upon constancy of /nl. In the following 
we adopted the large-sca le structure convention ( as opposed to 
the CMB convention. seelAfshordi & Tolleyll2008l: ICarbone et all 
l2008l|Pillepich. Porciani. & Hahnll2009l and lGrossi et alj|2009h "for 
defining the fundamental parameter /nl. According to this, the pri- 
mordial value of O has to be linearly extrapolated at z = 0, and 
as a consequence the constraints given on /nl by the CMB have 
to be raised b y ~ 30 per cent to comply with this pa per's conven- 
tion (see also IFedeli. M oscardini . & Matarre se 2009 for a concise 
explanation). 

In the case in which / N l + the potential O is a random 
field with a non-Gaussian probability distribution. Therefore, the 
field itself cannot be described by the power spectrum P<o(k) = 
Bk"~ 4 alone, rather higher-order moments are needed. The dom- 
inant higher-order contribution is generically given by the bispec- 
trum i?*(ki , k2, k 3 ). Only in very peculiar situations the bispectrum 
vanishes, and one has to resort to the trispectrum or higher-order 
correlations. The bispectrum is the Fourier transform of the three- 
point correlation function (0(ki)®(k 3 )0(k 3 )) and it can hence be 
implicitly defined as 

(cDtkOOKka)*^)) = (2/r) 3 <5 D (k 1 +k 2 +k 3 )S (ki,k 2 ,k 3 ). (2) 

As mentioned above understanding the shape of non- 
Gaussianity is of fundamental importance in order to pinpoint the 
physics of the early universe and the evolution of the inflaton field 
in particular. For this reason, in this work we considered four dif- 
ferent shapes of the potential bispectrum, arising from different 
modifications of the standard inflationary scenario. They are all de- 
scribed in the following. 



Local shape 

The standard single-field inflationary scenario generates negligibly 
small deviations from Gaussianity. These deviations are said to be 
of the local shape, and the related bispectrum of the Bardeen's po- 
tential is maximized for squeezed configurations, where one of the 
three wavevectors has much smaller magnitude than the other two. 
In this case the par ameter / N l must be a constant (see however 
iBvrnes et alj|2010"3 lbl). and it is expected to be of the same orde r 
of the slow-roll parameters dFalk. Rangaraian. & Srednicki|[l993h . 
that are very close to zero. 

However non-Gaussianities of the local shape can also be 
generated in the case in which an additional light scalar field, 
different from the inflaton, contributes to the observe d curvature 
perturbations dBabich. Creminelli. & Z aldarriaga 2004 1). This hap- 
pens, for instance, in curvaton models dSasaki. Valiviita. & Wands! 
120061 : lAssadullahi. Valiviita. & WandJ 120071) or in multi-fields 
mode l s dBartolo. Matarrese. & Riottol 12003 : iBernardeau & Uzanl 
2002; Huang 2001). In this case the parameter / NL is allowed to 
be substantially different from zero, and the bispectrum of the pri- 
mordial potential assumes the simple form 



Bo(ki, k 2 , k 3 ) = 2/nlB 2 [*r 4 C 4 + K"K A + ■ 



(3) 



O = Oq + /NL * (*G - <*G>) • 



(1) 



Due to its simplicity, the local model is the most studied one, 
especially in terms of cosmolo gical simulations (see however 
Wagner, Verde. & Bou bekeur 2 010) . 
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Equilateral shape 



Orthogonal shape 



In some inflationary models the kinetic term of the inflaton 
Lagrangian is not standard, containing higher-order derivatives 
of the field itself. One significant example of this is the DBI 
mode l dAlishahiha. Silverstein. & Tond 120041: Isilverstein & Tond 
2004 see also Arkani-Hamed et alj |2004l : ISeerv & Lidsevl 12005 ; 
Li. Wang. & Wand 120081) . In this case the primordial bispectrum 
is maximized for configurations where the three wavevectors 
have approximately th e same amplitude, and it takes the form 
dCreminelli et al.ll2007h 
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The function y(ki , k 2 , k 3 ) in the first line of the previous equa- 
tion represents a running of the parameter / NL that we have the lib- 
erty to insert since this parameter is not forced to be constant in 
the present case. This running has been considered as an actual part 
of the equil ateral bispectrum in all the calculations that fol low. It 
reads dChenll2005ljLo Verde et alj2008l ; ICrociani et alj|2009h 



7(k 1 ,k 2 ,k 3 ) : 



k\ + ki + £: 3 



(5) 



We adopted the exponent k = -0.2, that increase the level of non- 
Gaussianity at scales smaller than that corresponding to k C MB = 
0.086/; Mpc~'. This coincides with the larger multipole used in th e 
CMB analysis by the WMAP team dKomatsu et al]|2009l . l2010h , 
C ~ 700. When referring to the equilateral shape in the rest of this 
paper we always mean the bispectrum given by Eq. 10 with k = 
-0.2, unless otherwise noted. 



Enfolded shape 

For deviations from Gaussianity evaluated in the regular Bunch- 
Davies vacuum state, the primordial potential bispectrum is of 
local or equilateral shape, depending on whether or not higher- 
order derivatives play a significant role in the evolution of the in- 
flaton field. If the Bunch-Davies vacuum hypothesis is dropped, 
the resulting bispectrum is maximal for squashed configurations 

dChen et alfew^lHolman & Tolley||2008h. 

Meerburg, van der Schaar, & Corasaniti d2009l) found a tem- 



plate that describes very well the properties of this enfolded- shape 
bispectrum (see also lCreminelli et alfcoioh , that reads 



S (k, , k 2 , k 3 ) = 6/ NL B 2 [k"- A ^- A + k\-%- 4 + k'^ 4 k';- 4 + 

+ 3(*r 4 cX 4 ) 2/3 - 

/,( n -4)/3,2(n-4)/3p,-4 , , (n-4)/3 , 2(n-4)/3 t *-4 , 
I a.^ 2_ 3 3 \ 2. 

,( n -4)/3,2(n-4)/3„,-4 , (n-4)/3 , 2(n-4)/3„,-4 
i 2 3 1 l 3 2 

+ 4"- 4)/3 fcf- 4)/3 ^- 4 + 4- 4)/3 4 ( "- 4)/3 fcT 4 )] . (6) 

No running of the / NL parameter was introduced for both the en- 
folded and orthogonal shapes. The reasons are discussed further 
below. 



A shape of the bispectrum can be constructed that is 
nearly orthogonal to both the loc al and equilateral forms 
dSenatore. Smith. & Zaldarriagal 120100 . In this case the potential 
bispectrum is well approximated by the following template 



flo(ki,k 2 ,k 3 ) 



6/ NL S 2 [3(4"- 4)/3 4 < "- 4,/3 /f 3 '- 



+ k: 



Xn-4)/3,2( n -4)/3,„-4 , , (n-4)/3 , 2(n-4)/3 , „-4 



"-3 ""I "2 T " 2 ""3 
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"T" A. | A.^ /v ^ "T" *V j 3 

+ 4"- 4)/3 4 ( "- 4,/3 ^- 4 ) - 8 (y^- 4 ^- 4 ^- 4 ) 2 ' 3 - 

- 3 (/tr 4 fcT 4 + /fcf 4 &f 4 + k'^k'f 4 )] . (7) 

Constraints on the level of non-Gaussianity compatible with the 
CMB in the local, equilateral and orthogonal scenario s were 
recently given by the WMAP team (Kom atsu et alJ[20~LC?) , while 
constraints on enfo l ded n on-Gaussianity from LSS were given by 
IVerde & Matarresd d2009n 

Although there is no theoretical prescription against a 
running of the /nl parameter with the scale in the enfolded and or- 
thogonal shapes, we decided not to include one. The reason for this 
is that there is no first principle that can guide one in the choice of 
a particular kind of running, and until now no work has addressed 
the p r oblem of a running for these shapes dFergusson & Shellardl 
2009; Ferg usson. Liguori. & Shellardl 2010). As can be noted, in 
all non-local cases except the enfolded one we defined the level of 
non-Gaussianity by equating the bispectrum normalization in the 
equilateral li mit to the same quantity for the local shape (see the 
discussion in lFergusson & S hellard 200^). 



3 HALO MASS FUNCTION AND LINEAR BIAS 

Primordial non-Gaussianity produces modifications in the statistics 
of density peaks, resulting in differences in the mass function of 
cosmic objects and the bias of dark matter halos with respect to the 
underlying smooth density field. In the following we summarize 
how these modifications have been taken into account in the present 
work. 



3.1 Mass function 

For the non-Gaussian modification to the mass function of cos - 
mic objects we adopted the prescription of iLo Verde et alj d2008l) . 
The main assumption behind it is that the effect of primor- 
dial non-Gaussianity on the mass function is independent of the 
prescription adopted to describe the mass function itself. This 
means that, if n ( ^ s } (M, z) and nps(M, z) are the non-Gaussian and 
Gaussian mass functions , respectively, computed according to the 
IPress & SchechteJ d 19741) recipe, we can define a correction factor 
R(M,z) = n ( ^(M,z)/n PS (M,z). Then, the non-Gaussian mass func- 
tion computed according to an arbitrary prescription, n(M,z) can 
be related to its Gaussian counterpart through 



n(M,z) = ft(M,z)n {lj> (M,z) 



(8) 



In order to compute «ps(M, z), and hence K(M,z), 

ILo Verde et al.l 1 20081) per formed an Edgeworth expansion 
dBlinni kov & Moessnei] 1 19981) of the probability distribution for 
the smoothed density fluctuations field, truncating it at the linear 
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Figure 1. Left panel. The normalized skewness in units of /nl as a function of mass, for different bispectrum shapes as labeled. In the orthogonal shape case 
the skewness for /nl > is negative, thus the absolute value is plotted. Right panel. The scale-dependent part of the non-Gaussian correction to the linear bias 
in units of /ml, f° r a 10 l4 M Q /i~' halo mass. Different line types refer to different shapes of the primordial bispectrum, as labeled. Note that the correction for 
the orthogonal shape is negative (for positive /nl), thus we plotted the absolute value. As specified in the text, here and for all subsequent calculations the 
equilateral bispectrum shape comprehend the running of /ml- 



term in cr M - The resulting IPress & Schechteil dl974l) -like mass 
function reads 



exp 
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(9) 



In the previous equation p m = 3HQ£l m0 /%nG is the comoving mat- 
ter density in the Universe, <r M is the rms of density fluctuations 
smoothed on a scale corresponding to the mass M, and c5 c (z) = 
Ac/D + (z), where D + (z) is the growth factor and A c is the critical 
linear overdensity for collapse. The function S 3 (M) = /^(AfJ/cr^ 
is the reduced skewness of the non-Gaussian distribution, and the 
skewness j-i 3 (M) can be computed as 



M R (k 1 )M R (k 2 )M R (k 3 ) x 
dk l dk 2 dk 3 



x «6(k,)0(k 2 )©(k 3 )>- 



(2nf 



(10) 



The last thing that remains to be defined is the function M R (k), 
that relates the density fluctuations smoothed on some scale R to the 
respective peculiar potential, 

2 T(lc)k 2 

M R {k)=--^—W R (k), (11) 

n il m ,Q 

where T(k) is the matter transfer function and W R (k) is the Fourier 

transform of the top-hat window fun ction. 

In this work we adopted the iBardeen et ail Jl986l) matter 
transfer function, with the shape factor correction of Sugivamj 



( 1995) . This reproduces fa irly well the more sophisticated recipe 
of lEisenstein& Hul (l998) except for the presence of the baryon 
acoustic oscillation, that anyway is not of interest here. We ad- 
dit ionally adopted as refe rence mass function th e prescription 
of ISheth & Toi-menl 420021) (see Ijenkins et alj|200ll ; [Warren et al.l 
2006; iTinker et al . 2008 for alternative prescriptions). Other ap- 
proaches also exist for computing the non-Gaussian correc- 
tion to the mass function, that give results is broad agree- 
ment with those adop ted here jMatarrese. Verde. & Jimene3l2000l : 
lD'Amicoetalj2 010). These semi-analytic prescriptions are known 
to be in good agreement with cosmologi cal simulations having 
local non-Gaussian initial conditions dGrossi et alj 1 20071 120091 : 
iDesiacques & Se liak 2010), while more recently the same has been 
proven to be true also for more generic primordial bispectrum 
shapes (Wag ner. Verde. & Boubekeu r 2010). 

In the left panel of Figure[TJwe show the mass dependence of 
the reduced skewness for the four non-Gaussian shapes that have 
been detailed above. It is interesting to note that in the orthogonal 
case the skewness for positive / NL would be negative, thus bring- 
ing to a reduction in the abundance of massive objects. Also, the 
equilateral shape is the only case in which S3(M) is not mono- 
tonic. In computing the non-Gaussian corrections to the mass func- 
tion we have taken into account the modification to th e critical 
overdensity for collapse su ggested bv lGrossi et ai] j2009h (see also 
Maggi ore & RiottduOlOallbTlcj) . according to which A c -> A c yg, 
with q ~ 0.8. 



3.2 Halo bias 

Recently much attention has been devoted to the effect of 
primordial non-Gaussianit y on halo bias, a nd the use thereof 
for co nstraining / NL foalal et al] 120081; IVerde & Matarresd 



20091; ICar bone. Verde. & Matarresel [2008). In particular. 



Matar rese & Verdd 1 20081) have shown that primordial non 



Constraints on primordial non-Gaussianity from future surveys 5 



Gaussianity introduce a scale dependence on the large scale halo 
bias. This peculiar ity allows to place already stringent constra ints 
from existing data dSlosar etal . 2008; Afshordi & Tolle\l2008h . 

The non-Gaussian halo bias can be written in a relatively 
straightforward way in terms of its Gaussian counterpart as 
dCarbone. Mena. & VerdekoiCh 



b(M,z,k) = b iG \M,z)+Mk)cr 2 M [b iG) (M,z) - l] 



(12) 



where the function j3 R (k) encapsulates all the scale dependence of 
the non-Gaussian correction to the bias, and reads 



— — fM R (Ox 
,M R (k) Jo 

B m (£, y/a,k) 



dc. 



(13) 



where a = k 2 + t 2 + 2k(/u. In the simple case of local bispectrum 
shape it can be shown that the function f3 R (k) should scale as oc k~ 2 
at large scales, so that a substantial boost (if /nl > 0) in the halo 
bias is expected at those scales. 

In the right panel of Figure Q] we show the scale dependence 
of the function fi R {k) f or a fixed halo mass. In the r emainder of this 
paper we adopted the Shet h. Mo. & Tormenl ( feOOlh formula for the 
Gaussian bias. Thus, since we express the correction to the halo 
bias as a function of the Gaussian bias itself, we did not need to 
take into account t he ellipsoidal collapse correction suggested by 
iGrossi et al.1 d2009h . 



4 HALO MODEL OF THE LARGE-SCALE STRUCTURE 

The halo model dMa & Frvl2000l ; ISeliakl2000h is a physical frame- 
work that allows the description of the correlation function of dark 
matter as well of different tracers of the LSS such as galaxies ad 
galaxy clusters. It is based on the assumption that the power spec- 
trum of particles (either dark matter particles or tracers) is given 
by the sum of two contributions: particle pairs residing in the same 
structure, and particle pairs residing in two different structures. The 
implicit hypothesis underlying this assumption is that no particles 
are found outside bound structures. Accordingly, the cross spec- 
trum P xy (k, z) of two different kinds of particles x and y (if x = y 
then P xx (k, z) is the power spectrum of particle type x) can be writ- 
ten as the sum of two terms, 



P, y (k,z) 



Pxy,i(k, z) + P x ,, 2 (k,z). 



(14) 



The first one is named 1-halo term, while the second is the 2-halo 
term. From the discussion above, it is immediately obvious that 
the first term dominates on very small scales, while the second is 
dominant at large scales. 

The exact form of the two terms depend on the exact kind 
of particle that we are considering. However, a common feature is 
that the 2-halo term should depend upon the bias of dark matter 
halos, sinc e it should represents p article pairs residing in separated 
structures dCoorav & S heth 2002). 



4.1 Galaxy spectrum 

In this case we set x = y = g, and the terms contributing to the 
galaxy power spectrum can be written as 



PggAk,z) = 
and 



Jo 



(N g (N g - 1)|M> . ,„ 
n(M,z) \u g (M,z,k)\ P dM (15) 



n\(z) 



P gg , 2 (k,z) = P L (k,z)<pl(k,z), 
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Jr»+oo 
n(M, 
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z)b(M,z,k) 



(N g \M) 
» g (z) 



u„{M,z,k)dM. 



(16) 



(17) 



Note that we have included a scale dependence in the halo bias, 
in order to account for the effect of primordial non-Gaussianity. 
The functions (N g \M) and (N s (N g - 1)|M) are the first and second 
moment respectively of the halo occupation distribution p(N g , M), 
which represents the probability that N g galaxies reside inside a 
dark matter halo of mass M. The function P L (k, z) = Ak n D 2 + {z)T 2 {k) 
is the linear matter power spectrum. 

We observe that the previous integrals do not actually extend 
down to M = 0, since there is a minimum halo mass M g below 
which no galaxy formation is possible, i.e., p(N g , M) = for M < 
M g . Additionally, not all galaxies at all redshifts are accessible to 
observations, hence M 5 should be the minimum galaxy halo mass 
that, at a given redshift, enters in our fiducial catalogue. We come 
back on the issue of the minimum mass further below. The function 
n g (z) is the mean number density of galaxies, and reads 



»g(z) 



o 



(M,z)(N g \M)dM. 



(18) 



The quantity u g (M, z, k) in Eqs. l |15l > and dl7t represents instead the 
Fourier transform of the galaxy number density inside dark matter 
halos of mass M and redshift z- We set it equal to the number den- 
sity of dark matter particles, that is u B ( M,z,k) = p(M,z,k)/M. For 
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Navarro. Fre nk. & White 



the da rk matter density profile we adopt a 

( 1996) shape (NFW henceforth, see also 

1 19951 1 19971) with a concentration- mass relation that reproduces 
simulated matter power sp ectra (see iHuffenberger & Selialdlixii 
iFedeli & Moscardinil2010] for discussions). We discuss more on the 
adopted density profile in Section|6] 

We make the assumption that if a dark matter halo hosts 
at least one galaxy, then one galaxy should sit at the center of 
the halo itself. As a consequence, the 1-halo term of the galaxy 
power spectrum consists of two contributions: galaxy pairs that in- 
volve the centra l objec ts and all the other pairs. As discussed by 
ICoorav&Sh"ethl (2002) we can self-consistently assign p = 1 to 
the first class ((N g (N s - l)\M) < 1) and p = 2 to the second 
({N g (N g - l)\M) > 1). As for the moments of the halo occupa- 
tion distribution we adopted simple fitting forms that reproduce the 
results of semi-analytic galaxy formation models, that is 



M 

{N g \M) = N g , \ — 



and 

{N g {N g - 1)|M> = f 2 (M)(N s \M) 2 



(19) 



(20) 



dCoorav & S he th 2002; Cooravl2004rHutsi & Lahavl2008l) . where 
/(M) = log VM/ICMq/i- 1 if M ^ 10 l3 M o /r' and /(M) = 1 
otherwise. This corresponds at assuming a Poisson distribution at 
high halo masses with deviations thereof occurring at low masses. 

The three free parameters N gl o, Mp and 9 depend on th e type 
of galaxy that is under consideration dCoorav &~S heth 2002). Typ- 
ically, for blue galaxies yV„ = 0.7, M = 4 x 10 12 M Q /r', = for 
M ^ Mo and = 0.8 otherwise. For red galaxies we have instead 
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Figure 2. The power spectra of galaxies and galaxy clusters obtained ac- 
cording to our halo model at z = 0, compared to the matter power spectrum 
and to the cluster-galaxy cross spectrum at the same redshift, as labeled. 

N g , a = 1, M = 2.5 x 10 12 M o /i-\ and 6 = 0.9. Obviously, the av- 
erage number of galaxies inside a halo of some mass M is the sum 
of the average number of red galaxies and of blue galaxies. Now, 
the distinction between blue and red galaxies can be rather quali- 
tative, and dependent on the framework we are interested in (e.g. 
observations vs. semi-analytic models). The slitless spectroscopic 
instrument planned for Euclid will be mainly sensitive to the flux 
and equivalent width of emission lines, mainly Ha, and thus will 
mostly observe star-forming galaxies. Although this category in- 
cludes dust-obscured starburst objects, that are not properly blue, 
we believe the parameter set for blue galaxies to be more suited 
for the target of Euclid. Hence we shall adopt it in the remainder 
of this work. In Section [6] we discuss how our results are changed 
if we consider a multi-slit spectroscopic instrument, that select a 
different galaxy population. 

4.2 Cluster spectrum 

For galaxy clusters the situation is less elaborated, since it is com- 
monly assumed that only one cluster is contained inside each dark 
matter halo that is massive enough (this is obviously not true for 
galaxies and dark matter particles). Consequently, the 1-halo term 
vanishes, P cc ,i(k, z) = 0, while the 2-halo term reads 



P cc , 2 (k,z) = P L (k,z)<pt(k,z), 



where 



tfc(k, z) 



Jn 
o 



b(M,z,k) 
(M,z) - dM. 



(21) 



(22) 



JO «cfe) 

In this case as well there is a minimum halo mass M c below 
which no cluster is formed (or, below which clusters are not ob- 
servable), and the integrals are to be evaluated above that mass. 
The function n c (z) then reads 



«c(z) : 



Jn 
o 



(M, z)dM. 



(23) 
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Figure 3. The minimum mass as a function of redshift for spectroscopi- 
cally selected galaxies (solid green line) and weak lensing selected clus- 
ters (dashed red line) adopted in the present work. The points are the av- 
erage galaxy halo masses derived from semi-analytic evolutionary models 
for two different Ha flux thresholds, log(fu a ) = -15.4 (red triangles) and 
\og(ffi a ) = -15.5 (blue circles, CGS units). 



Further below we detail how the issue of the minimum cluster mass 
is addressed in this paper. 

In Figure [2] we show the power spectrum of galaxy clus- 
ters alongside that for galaxies, and compare them to the 
dark matter power spectr um obtained with the halo model 
dFedeli & Moscardinfo lO). At linear scales clusters are more bi- 
ased than galaxies as is to be expected. Moreover, since the power 
spectrum of clusters is made only by the 2-halo term, their corre- 
lation drops off substantially at small scales as compared to other 
tracers. 



4.3 Cluster-galaxy cross spectrum 

Extending the previous results. iHiitsi & Lahav (2008) derived the 
two halo model contributions to the cross correlation of clusters 
and galaxies. They can be written, respectively, as 



Pc g ,i(k,z) 



and 



Jo 



(NAM) 

n(M, z) ' uJM, z, k)dM 

n g (z)n c (z) " 



P, g 2(k,z) = PL(k,z)<Pc(Kz)<p g (k,z). 



(24) 



(25) 



In Figure|2]we also show the cross spectrum of galaxy clusters 
and galaxies. As can be seen, at large scales the cross spectrum 
stays in between the spectra of galaxies and clusters. However, it 
shows more power on small scales than both clusters and galaxies. 
The fact that the scale dependence of the cross spectrum is different 
from that of the galaxy spectrum and cluster spectrum means that 
it carries a different kind of information on the spatial distribution 
of LSS tracers, thus the cross-correlation can make up an important 
improvement when constraining cosmological parameters. 
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4.4 Minimum masses 

As mentioned above, in order to compute the cluster and galaxy 
power spectra it is necessary to select a minimum mass for both 
classes of objects. In general thi s minimum mass depends on red- 
shift. Let us start with clusters. IBerge. Amara. & Refregier d2010h 
have estimated the weak lensing selection function for galaxy clus- 
ters given a Eu clid-like survey and different S/N detection thresh- 
olds (see also [Laureiis 2003). In order to be conservative, we 
ad opted S/N ^ 5, and fitted t h e corr esponding curve in Figure 1 
of IBerge. Amara. & Refregierl 120101) . The result is displayed by 
the red dashed curve in Figure [3] which represents the M c (z) we 
adopted in this work. 

For galaxies, we adopted the minimum halo mass M g (z) = 
10"M Q /i~' which, given the halo occupation moments reported in 
Eqs. l |19t and l |20| >, reproduces in the reference Gaussian cosmol- 
ogy the effective bias predicted for Ha selected galaxies by Euclid 
dOrsi et al] |2010). According to Eq. J 1 6b . the effective galaxy bias 
can be estimated within the framework of the halo model as the 
large scale limit of tp g (k, z), i.e. b e (z) = lim^ p (p g (k, z). The effective 
bias estimates reported in lOrsi et al.l (|20J_Oj) have bee n performed 
adopting semi-analitic m odels of galaxy formation teaugh et al.l 
l2005tlBoweretalJl2006h . and assuming a spectroscopic selection 
based on Ho- emission. The latter is appropriate for our scope, since 
we are considering star-forming galaxies o nly. The points re ported 
in Figure[4] are the pred ictions performed bv lOrsi et alj bold) using 
the lBaughetai] ( l2005h model for two different Ha flux thresholds 
(log(/ Htt ) = -15.4 and log(/ H „) = -15.5 in CGS units). Obviously 
there is not much difference between the two choices, that are well 
reproduced by the solid line, representing the outcome of our halo 
model. In Figure[3]we also report the minimum mass of galaxy ha- 
los, co mpared to the average such mass found again bv lOrsi et al.l 

feoich . 



5 RESULTS 

In this Section we summarize our main results. First, we explored 
what happens to the power and cross spectra described above when 
the distribution of primordial density fluctuations deviates from a 
Gaussian. Then, we performed a Fisher matrix analysis in order to 
forecast the constraints on the main cosmological parameters that 
can be expected by a Euclid-like all-sky survey. 

5.1 Effect of non-Gaussianity on the power and cross spectra 

As a first step of our analysis we computed the Euclid-expected 
galaxy power spectrum, cluster power spectrum, and cluster-galaxy 
cross spectrum for all the cosmological models considered in this 
work, and investigated the modifications induced by primordial 
non-Gaussianity on each one. In Figure [5] we show the effect of 
the different shapes of primordial non-Gaussianity examined in the 
present paper on the power and cross spectra for different tracers of 
the LSS. We show results for / NL = ± 100, since this roughly corre- 
sponds to the best high confidence limits that are currently placed 
on this parameter, and for /nl = ±200 as a comparison. We recall 
that the power spectrum of dark matter halos in cosmological sim- 
ulatio ns with local non-Gaussian initial conditions has been stud - 
ied in lGrossi et all J2009l) and lDesiacques. Seliak. & Ilievl J2009t) . 
However, here we are interested in a wider variety of bispectrum 
shapes, and in specific classes of objects that are not necessarily in 
a one-to-one relation with dark matter halos. 




i i i i I i i i i I i i i i I i i i i I i i i i 
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Figure 4. The effective bias as a function of redshift for spectroscopically 
selected galaxies (solid black line). The points are the galaxy biases derived 
from semi-analytic evolutionary models for two different Ha flux thresh- 
olds, log(fu a ) = -15.4 (red triangles) and log(fn a ) = -15.5 (blue circles, 
CGS units). 



The corrections to the power and cross spectra that are due to 
different kinds of primordial non-Gaussianity resemble the correc- 
tions to the large scale halo bias that have been shown in Figure[T] 
This is to be expected, since in the halo model the mass function is 
normalized over (see Section^, hence the largest effect is expected 
to be due to the halo bias. In particular, for the local shape we ob- 
serve a sharp increase of the power at large scales for all tracers. 
A similar increase in power at large scales, although not as strong, 
is also observed for the enfolded and orthogonal shapes. Note also 
that the trend for the latter is reversed, due to the fact that in this 
case a positive /nl implies a negative effect on the bias and the 
mass function. Finally, for the equilateral shape we find a slight in- 
crease in the power at intermediate scales, again in agreement with 
the behavior of the halo bias. 

As one might expect, the larger effect is found in correspon- 
dence of the tracers that are more biased with respect to the underly- 
ing matter density field, that are galaxy clusters. This can be under- 
stood by looking at Eq. &121 , where it is shown that the difference 
between the non-Gaussian halo bias and its Gaussian counterpart 
is in fact proportional to the square of the Gaussian bias itself. On 
the other hand, the effect on the cluster-galaxy cross power spec- 
trum is at an intermediate level between the effect on the galaxy 
power spectrum and that on the cluster power spectrum. Due to 
their low average bias, star-forming galaxies are very poorly af- 
fected by primordial non-Gaussianity, and this is particularly true 
for the enfolded and orthogonal bispectrum shapes. 

A closer look to the first and third column panels of Figure [5] 
reveals an interesting fact. While the change in the cluster-galaxy 
cross spectrum due to primordial non-Gaussianity is larger than 
the change in the galaxy power spectrum, it is more so for the en- 
folded and orthogonal bispectrum shapes, as compared to the local 
and equilateral ones. This means that the inclusion of the cross- 
correlation between clusters and galaxies should bring a better im- 
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Figure 5. The ratio of the power and cross spectra for different tracers of the LSS computed in various non-Gaussian cosmologies to the same quantities 
evaluated in the reference Gaussian case. From top to bottom: non-Gaussianities with local, equilateral, enfolded, and orthogonal bispectrum shapes. From 
left to right: the power spectrum of star-forming galaxies, galaxy clusters, and the cross spectrum of clusters and galaxies. For all non-Gaussian cosmologies 
considered here values of /nl = 100 (solid red lines), /nl = 200 (solid green lines), /nl = _ 100 (dashed red lines), and /nl = -200 (dashed green lines) 
were chosen. All results are shown at z = 0. 



provement over the constraints on non-Gaussianity obtained by us- 
ing the galaxy spectrum alone in the first two cases than in the 
second ones. Although it is difficult to gauge what the effect of the 
additional inclusion of the cluster power spectrum would be, we 
show in the next subsection that this expectation is to some extent 
satisfied. 



5.2 Fisher matrix analysis 

For the Fisher information matrix of the cross correlation be- 
tween clusters and galaxies we followed the calculations of 
iHtitsi & Lah av (2008), assuming perfect redshift knowledge for 
both galaxies and clusters. This assumption is justified since galax- 
ies are spectroscopically selected, while clusters, being chosen as 
the highest S/N cosmic shear peaks, should be very massive, thus 
making follow-up confirmations relatively straightforward. In ad- 



dition, we assumed redshift space distortions to be sufficiently well 
understood to be modeled away. Accordingly, the Fisher matrix for 
the redshift bin centered at z can be well approximated as 



Viz) 
4tt 2 



■ . „ d\nP cg (k,z)d\nP cg (k,z) 
k dkw cg {k,z)- 



where w cg Qc, z) = 2/ [l + y cg (k, and 



y cg (k, z) 



[1 +n c (z)P cc (k,z)][l +n,(z)P,Jk,z)] 
n c (z)n g (z)Pl g (k,z) 



(26) 



(27) 



In the above set of equations, the quantities n c (z) and n g (z) 
are the average number densities of clusters and galaxies in the 
survey at hand (Eqs.[23]and[T8] respectively), V(z) is the comoving 
volume contained in the unit redshift around z, and is the «-th 
parameter of our cosmological model. The two wavenumbers k m i n 
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Table 1. The adopted parameter set for our fiducial Euclid-like survey 



GALAXIES 
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Sky Coverage 

Redshift Efficiency 

Ha Flux Threshold 

Cosmic Shear S/N Threshold 
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Figure 6. The joint 68% confidence level contours in the /nl - erg plane given by different combinations of the galaxy and cluster power spectra and the 
cluster-galaxy cross spectrum, as labeled. Results are shown for all the four shapes of the primordial bispectrum that have been considered in this work. All 
other cosmological parameters have been kept fixed to their WMAP-7 fiducial values. 



and A; max represent the boundaries of the wavenumber range used in 
the analysis. The chosen redshift bins should be relatively narrow, 
in order to treat the spatial number densities of tracers within each 
bin as constants.The total Fisher matrix is then given by the sum 
over the n z adopted redshift bins, 



(28) 



In the case in which we are interested only in galaxies or only 
in clusters, the formalism remains the same, with the only differ- 
ence being given by the replacements of P cg (k,z) with P xx (k,z) 
(with x = c or x = g) and of w C6 (rc, z) with 



Wxx(£,Z) : 



n x (z)P x (k,z) 



l+n^z)P x (k,z) 



(29) 



in Eq. l !26t . In our analysis we adopted a redshift range typical 
for Ho- galaxies selected by Euclid, 0.5 ^ z ^ 2.1 (cf. the Eu- 
clid Yellow Book), and a range of scales corresponding to & m i n ^ 
k ^ 0.3 Mpc -1 , where k m i n matches the largest scale available for 
a given redshift bin. For completeness, in Table [TJ we summarize 
all the survey parameters that have been adopted in our Fisher ma- 
trix analysis including the redshift efficiency, that is the fraction 
of spectroscopically selected galaxies for which a reliable redshift 
measurement is expected. 

In Figure [6] we show the joint 68% confidence level contours 
in the / NL - <r 8 plane that result from the Fisher matrix analysis per- 
formed for all the four bispectrum shapes considered in the present 
paper. All other cosmological parameters are have been kept fixed 
to the fiducial WMAP-7 values summarized in Section [TJ We re- 
port the results of considering the galaxy power spectrum alone, 
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Table 2. The 1 - cr forecasted errors for /nl and erg . 
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the cluster power spectrum, the cluster-galaxy cross spectrum, and 
different combinations of the three. For further clarity, in Table [2] 
we report the 1 - cr errors that are forecasted both on the level 
of non-Gaussianity / NL and on the amplitude of the matter power 
spectrum erg . The fiducial model that we adopted has /nl = and 
<r 8 = 0.809, as specified in SectionQ] 

As can be seen, the tracers of the LSS that individually give 
the best constraints on both parameters considered are galaxies, due 
to their very large number density which beats down the shot noise 
in the Fisher matrix. Galaxies alone give constraints on erg that are 
more than an order of magnitude better than those of clusters alone, 
while the performance on the level of primordial non-Gaussianity 
/nl, as expected, depends on the specific choice of the primordial 
bispectrum shape. Galaxies perform a factor of ~ 8 better than clus- 
ters in the equilateral case, but only a factor of ~ 3 better in the 
orthogonal configuration. 

The constraints on the /nl - erg plane that we gather 
fr om clusters al o ne are somewhat different than those reported 
in ISartoris et al. which performed a similar Fisher ma- 

trix analysis on the cluster sample expected to be detected by 
the planned Wide Field X-r ay Telescope (WFXT) (see also 
ICunha. Huterer. & Dorel l20ld for a related work). Particularly, 
while the forecasted error on /nl is comparable, the error on cr 8 
is much smaller in our case than in their. We ascribe this differ- 
ence to two factors. Fir st, clusters are selecte d in different ways in 
the two cases. Second, ISartoris et alj d2010t) introduced nuisance 
parameters for estimating the uncertainty in the mass calibration, 
and at the same time adopted priors on several cosmological pa- 
rameters inspired by Planck forecasts. We infer the second factor 
to be t he dominant one, as rerunning the analysis of ISartoris et al.l 
d2010h without nuisance parameters and priors brings to errors on 
erg comparable with ours own (B. Sartoris, private communication). 
We believe it to be safe to assume that Euclid weak lensing clus- 
ter catalogs will have sufficient control over systematics to ignore 
the contribution of imperfect knowledg e of the mass-observable re- 
lation. This point of view is shared by Ber ge. Amara. & Refregier 

feoich . 

The behavior of the confidence contours given by the cluster- 
galaxy cross spectrum alone is interesting. Firstly, while the con- 
straining power on /nl is greatly enhanced as compared to the 
constraining power of clusters alone, the same is not true for <x 8 . 
Actually the I - cr errors on <x 8 expected from the cluster-galaxy 
cross correlation are systematically ~ 15% larger than those ex- 
pected from clusters alone. As a consequence, while galaxies still 
play a dominant role about constraining cr 8 , the cross correlation 
gives comparable constraints on /nl- As a matter of fact, for the 
enfolded and the orthogonal shapes, the cross spectrum performs 
even better with respect to /nl, so that combination of all the three 
probes brings to a significant improvement over galaxies alone. 



Specifically, while the combined probes decrease the 1 - cr error 
on /nl only by ~ 30% for the local bispectrum shape, this decre- 
ment reaches a factor of ~ 2 in the orthogonal case. The fact that 
the improvement in constraining power due to the inclusion of the 
cross spectrum is more pronounced for the enfolded and orthogo- 
nal shapes can be deduced from the observation that we made at the 
end of the last Subsection, justified by the complicated interplay of 
galaxy and cluster biases modified by primordial non-Gaussianity. 

A perhaps unexpected result is that the constraints on the / N l 
parameter in the equilateral case are almost comparable with those 
for the local case, and much more stringent than for the enfolded 
and orthogonal shapes. Since the effect of non-Gaussianity on the 
large-scale bias is very slight in the equilateral case one would 
have naively expected the opposite behavior. However, it should 
be recalled that the correction to the mass function (which also 
enters the calculations of the power spectra and the Fisher matri- 
ces) for the equilateral shape is comparable to that for the local 
shape. Moreover, in the former case the non-Gaussian correction 
to the bias is concentrated at intermediate scales, and the resulting 
effect on the power spectra is thus less degenerate with erg with 
respect to the other cases. Finally, the integral over wavenumbers 
in the definition of the Fisher matrix weights more smaller scales. 
Hence a modification of the power spectrum at a given level has a 
stronger impact if applied to sma l l scal es rather than large ones. In 
iFedeli, Moscardini. & Matarresd d2009h it is shown how the effect 
of equilateral primordial non-Gaussianity on the mass function and 
halo bias is reduced when removing the running of / N l- In the next 
Section we show that this has the effect of loosening the constraints 
on the level of primordial non-Gaussianity. 

As a final step of our analysis, we computed the correlation 
coefficients between /nl and erg for different non-Gaussian shapes 
and for different combinations of LSS probes adopted. The corre- 
lation coefficient of the j— th and the j—th cosmological parameters 
in our set is defined as 

nj = p== (30) 
V(/ : > (/ : ) 

where T is the Fisher matrix relative to the probe we are consid- 
ering. For instance, the Fisher matrix for the combination of the 
cluster power spectrum with the cluster-galaxy cross spectrum is 
T = ^cc + ^cg- According to its definition, the correlation coef- 
ficient would have value ry = +1 for perfectly correlated param- 
eters, Tjj = -1 for perfectly anti-correlated ones, and = for 
uncorrelated parameters. The results are summarized in Table [3] 
The first thing to note is that, while the coefficients are positive for 
the local, equilateral, and enfolded shapes, they are negative for the 
orthogonal shape. This reflects the point that, in the latter case, the 
confidence contours plotted in Figure [6] are tilted in the opposite 
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Table 3. The /nl-0"8 correlation coefficients. 



LOCAL EQUILATERAL ENFOLDED ORTHOGONAL 



P gg 0.05 0.02 0.05 -0.09 

Pec 0.37 0.33 0.27 -0.18 

P cg 0.13 0.13 0.08 -0.04 

P cc + P cg 0.13 0.12 0.08 -0.05 

Pgg + Pec + P cg 0.04 0.02 0.03 -0.04 



direction with respect to the former. This is a consequence of the 
fact that in the orthogonal case, a positive value of / NL corresponds 
to a negative skewness of the density field and a negative correc- 
tion to the linear bias, while for the other shapes considered in this 
work the opposite is true. A closer look to the absolute value of the 
correlation coefficients reveals that the constraints given by galaxy 
clusters alone are the most degenerate between /nl and as, while 
the constraints given by galaxies alone are almost undegenerate, in 
the sense that they pin down <x 8 much better than /nl- 

The positive correlation between <r g and / NL in terms of power 
spectrum of cosmic tracers can be understood by considering the 
effect that both parameters have on the halo bias. An increase in <r 8 
has the effect of making structures less rare, and hence less biased 
with respect to the background matter density field. On the other 
hand an increase in / NL in the positive direction increases the bias, 
except for the peculiar orthogonal case, that indeed shows an oppo- 
site degeneracy. The resulting degeneracy is however not very large 
because an increase in cr 8 also increases the matter power spectrum, 
thus partially counteracting the decrement in the halo bias. 



20 
10 



-10 
-20 



68% C.L. jointly 
f NL -equilateral, y= 
slitless 

S/N>5 



1 

i'V 


\ i 


A ■ 


. r. 

j j i 




;'/ 






K 

■ v.*. 





galaxy 

cluster 

cross 

dust er+ cross 



0.76 0.78 



0.80 

^8 



0.82 0.84 



Figure 7. The joint 68% confidence level contours in the /nl - erg plane 
given by different combinations of the galaxy and cluster power spectra and 
cluster-galaxy cross spectrum. We show only the results for the equilateral 
shape of the primordial bispectrum, assuming no running of the parameter 
/nl- 



6 DISCUSSION 

In this Section we debate some of the assumptions that have been 
adopted in the present paper, an d discuss overlap pings with other 
works. We begin by noting that Wan g et all J20ld) forecasted the 
constraints on the redshift evolution of the dark energy equation 
of state parameter expected by the galaxy power spectrum mea- 
sured by Euclid. The error on tr 8 that they found is of the order 
of ~ 7 x 1CT 2 , thus being much larger tha n our result . The rea- 
sons for this mismatch are the following: tl lWang etalj feOlCh let 
all the cosmological parameters free to vary, while we fixed them 
all to their WMAP best fits, except for /nl and cr 8 ; i i) we used the 
full power spectrum information, while lWang et alj d20ld) limited 
themselves only to the shape and position of the Baryon Acoustic 

Oscillation. 

Recentlv. lSmith. Desiacques. & M arian (2010) have shown by 
using numerical simulations of structure formation in the presence 
of local non-Gaussianity, that dark matter halos tend to be on aver- 
age more centrally concentrated than their counterparts in the ref- 
erence Gaussian cosmology for positive f NL , while the opposite is 
true f or negative /nl. This fin ding is in agreement with previous 
work jAvila-Reese et all2003h . and with the naive expectation that 
a higher efficiency in forming high-mass objects implies. Under 
our assumption that the galaxy number density traces the dark mat- 
ter density, this would underpin a different distribution of galaxies 
within dark matter halos as well. However, such an assumption is 
probably rough, and it is not clear if the galaxy distribution would 
be modified similarly to the matter distribution. Nevertheless, ef- 
fects due to the inner structures of clusters show up only on very 
non-linear scales, that are important for cosmic shear studies, but 



have been neglected in the present investigation (see Section [5]l- 
Thus, we stick to our choice of a standard NFW profile for the dis- 
tribution of galaxies within halos for all cosmologies considered 
here. 

The running of the /nl parameter introduced for the equilat- 
eral shape in Eq. © also deserves further discussion. Namely, dif- 
ferent authors tend to use different sha pes for this running, or, i n 
some circumstances, no running at all dVerde & Matarresa 20091) . 
Therefore, we redid our Fisher matrix analysis by setting the run- 
ning y(ki,k 2 ,k 3 ) = 1. The results are shown in Figure|7] As can 
be seen, while the constraints on cr 8 are almost unchanged those 
on /nl are visibly loosened. For instance, the 1 — er forecasted 
error obtained by galaxy clusters alone goes from ~ 13 to ~ 30 
when the running is removed, while that obtained by galaxies alone 
grows from ~ 1.5 to ~ 4.6. This kind of effect was expected, 
since the running we introduced had the effect of increasing the 
level of non-Gaussianity on sub-CMB (cluster) scales. By setting 
7(k!,k2,k 3 ) = 1 such an increase is removed, and the errors in- 
crease. This result stresses the importance of considering the same 
/nl running when results from differe nt works are compared (see 
also lShandera. Dalai. & Hutererll2010h . 

In order to make our analysis more complete, we evaluated 
constraints on non-Gaussianity considering survey configurations 
different from the fiducial one described in the preceding Sections. 
As a first step, we changed the S/N threshold for the detection of 
dark matter halos in cosmic shear maps, going from S/N ^ 5 to 
S/N ^ 7. We again fitted the cor responding curve in Figure 1 of 
Berge rAmara. & Refregiej j2010h in order to infer the minimum 
cluster mass to be adopted. In this case, M c (z) at z = is a fac- 
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Figure 8. The joint 68% confidence level contours in the /nl - erg plane given by different combinations of the galaxy and cluster power spectra and the 
cluster-galaxy cross spectrum, as labeled. Results are shown for the local (left panel) and orthogonal (right panel) bispectrum shapes, assuming a weak lensing 
cluster detection threshold of S/N ^ 7. 



tor of ~ 2 larger than for the S/N ^ 5 case. Only two bispec- 
trum shapes have been considered here: the local one, since it is 
the one expected to give the best constraints, and the orthogonal 
one, since it is the less studied one. In Figure [8] we show the joint 
68% confidence levels on <r 8 and / NL given by this new detection 
criterion for galaxy clusters, and in Table [4] we report the numeri- 
cal values of the respective errors. While the constraints given by 
galaxies alone are obviously unchanged, those given by clusters 
and the cluster-galaxy cross-correlation are significantly loosened. 
As a consequence, very little information is added if combining the 
galaxy power spectrum with power spectra involving clusters, and 
this remains true for the orthogonal model as well. This result im- 
plies that, while selecting cosmic shear peaks with higher S/N we 
are effectively selecting more massive structures, on whose statis- 
tics the effect of non-Gaussianity is larger, the respective constrain- 
ing power is much too suppressed by the exponential decrement in 
the abundance of such objects. 

As a second alternative survey configuration we considered 
a multi-slit spectroscopic instrument on the model of the Digital 
Micro-mirror Devic es (DMDs) conceived for the SPACE mission 
dCimatti et ai1l2009T) . This is particularly relevant in view of future 
space missions alt ernative to Euclid, such as WFIRST and JEDI 
dCrotts et alj|2005t) . In this case selection of galaxies would occur 
according to their H-band flux. As a consequence Luminous Red 
Galaxies (LRGs) would be detected as well, thus we had to use the 
halo occupation distribution relative to a mixture of red and blue 
galaxies, rather than to blue galaxies only as we did for our fidu- 
cial survey. We found that a new constant minimum halo mass of 
M g (z) = 3 x 10"Afoi"' reproduces fairly well the effective bias of 
galaxies with H-band m agnitude br i ghter than H AB = 22 accord- 
ing to the simulations of lOrsi et al.l d201fj) . All the other parame- 
ters of the survey are left unchanged, while the redshift range of 
the analysis is extended to 0.1 ^ z ^ 2.1 (for both clusters and 
galaxies, although the former are unaffected by spectroscopy) in 
order to comply with the different observational specifications. In 
Table [4] we report the 1 - <x errors forecasted for this alternative 
configuration. As can be seen, in this case the constraints on both 
/nl and erg are significantly improved with respect to our fiducial 
Euclid-like case. The error on the level of non-Gaussianity can for 
instance reach A/ NL ~ 0.5 for the local bispectrum shape. Since 
the cluster selection is unchanged, we interpret this improvement 



as due to a combination of the different galaxy selection and of the 
extended redshift range that is allowed by the multi-slit configura- 
tion. We expand a little bit more on this further below. This im- 
provement is additionally emphasized in Figure [9] where we show 
the 1-parameter 68% confidence levels on the <x 8 -/ NL plane given 
by the combination of the power spectra of clusters and galaxies 
and the cluster-galaxy cross spectrum for both the multi-slit con- 
figuration we are considering and the fiducial slitless Euclid case. 
It is interesting to note that, due to the different galaxy population 
selected in the multi-slit case, the inclination of the confidence el- 
lipses changes with respect to the fiducial slitless configuration. 

Finally, we estimated how the constraints on the level of non- 
Gaussianity changed upon modification of the redshift and scale 
ranges considered in the Fisher matrix analysis, for both the slitless 
Euclid-like configuration and the multi-slit case. In Figure [TT)1 we 
show the results of changing either the minimum scale included in 
the analysis, k m . m , the minimum, and the maximum redshift of both 
galaxies and clusters. Fiducial vales are k miix = 0.3 Mpc -1 , 0.5 
z ^ 2. 1 for Euclid, and 0.1 ^ z ^ 2.1 for the multislit case. As 
one could naively expect, the error on / NL decreases by increasing 
both the minimum scale and the maximum redshift (although in the 
latter case the trend is quite mild), while it increases by increasing 
the minimum redshift (because the maximum redshift is held fixed). 
As mentioned above, the differences between the slitless and the 
multi-slit case are due to a combination of the different selection 
of galaxies with the different redshift range adopted. An exception 
to this is given by the two middle panels of Figure [10] where the 
redshift ranges adopted are the same for both configurations, and 
thus the differences can be ascribed only to the different selection 
functions. 



7 SUMMARY AND CONCLUSIONS 

We adopted the physically motivated halo model in order to com- 
pute the effect of different kinds of primordial non-Gaussianity 
on the power spectrum of galaxies and galaxy clusters, as well as 
on the cluster-galaxy cross spectrum. Specifically, we considered 
galaxies selected spectroscopically according to their Ho- flux, and 
galaxy clusters selected as the largest S/N peaks in cosmic shear 
maps, having in mind future wide field optical/near-infrared sur- 
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Table 4. The 1 — cr forecasted errors for /nl and erg assuming alternative survey configurations. 
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Figure 9. The 1 -parameter 68% confidence level contours in the /nl - cr 8 plane given by the combination of the galaxy power spectrum, the cluster power 
spectrum, and the cluster-galaxy cross spectrum. Results are shown for the local (left panel) and orthogonal (right panel) bispectrum shapes. We illustrate the 
comparison between the outcome for our fiducial Euclid-\ike survey and for the multi-slit case discussed in the text, as labeled. 



veys such as Euclid and WFIRST. We additionally performed a 
Fisher matrix analysis in order to forecast the errors on the joint 
determination of the level of non-Gaussianity /nl and the ampli- 
tude of the matter power spectrum er g . The main findings of this 
work can be summarized as follows. 



• The non-Gaussian corrections to the power spectrum of tracers 
of the LSS resembles the modifications to the large-scale bias of 
dark matter halos, as one might naively expect. The largest effect is 
seen for the local bispectrum shape, while the smallest appears for 
the enfolded shape. The power spectrum of massive galaxy clusters 
is more heavily modified than the spectrum of galaxies, because the 
former are originally more biased than the latter. The modification 
to the cluster-galaxy cross spectrum lies somewhere in between. 

• Galaxies have a much higher constraining power on both /nl 
and erg as compared to galaxy clusters, due to the much lower abun- 
dance of the latter that is not adequately compensated by the larger 
effect on the relative power spectrum. Assuming a Euclid-like sur- 
vey, while spectroscopically selected galaxies can constrain /nl to 
the level of a few and erg to the level of ~ 3 x 1CT 4 , errors on pa- 
rameters derived by clusters alone are at the level of A/nl ~ 10 and 
Acrg ~ 8 x 10~ 3 (with some dependence on the primordial bispec- 
trum shape). 

• When considering the cluster-galaxy cross spectrum alone, the 
forecasted constraints on erg are comparable to the cluster-only con- 
straints, while the constraints on tr 8 are improved by more than one 
order of magnitude, reaching a predicted error comparable with that 
coming from galaxies alone, A/nl ~ a few. This result highlights 



the high complementarity of the cluster power spectrum and the 
cluster-galaxy cross spectrum as cosmological probes. 

• While the constraints on erg coming from the galaxy power 
spectrum alone are almost unchanged when combined with the 
cluster-galaxy cross correlation, the constraints on / NL are im- 
proved. This is true only to a slight level for the local and equi- 
lateral bispectrum shapes, however the error on / N l can be reduced 
by a factor of ~ 2 for the enfolded and orthogonal cases. The addi- 
tion of the cluster power spectrum carries only slight change to this 
conclusion. 

• As expected the parameters / NL and er 8 are degenerate with re- 
spect to the power spectra of LSS tracers. This degeneracy is more 
marked for the cluster power spectrum, with correlation coefficients 
reaching up to 0.3 - 0.4. The degeneracy on the other hand is al- 
most insignificant for the galaxy power spectrum, in the sense that 
the constraining power weights much more on tr 8 than on / N l- The 
degeneracy with respect to the cluster-galaxy cross spectrum lies 
in between the former two, with the exception of the orthogonal 
bispectrum shape. 

• We considered several survey configurations alternative to the 
fiducial Euclid-like one, finding that a multi-slit spectroscopic in- 
strument would allow more stringent constraints both on / N l and 
er 8 . This improvement is due to a combination of the fact that with 
this other arrangement the selected galaxy population is different, 
and that the data analysis can be pushed down to a lower minimum 
redshift. This result is interesting with respect to future wide field 
survey concepts alternative to Euclid, such as WFIRST. 
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Figure 10. The 1 - a forecasted error on the level of non-Gaussianity given by the combination of the cluster and galaxy power spectra and the cluster-galaxy 
cross spectrum for both the Euclid and the multi-slit survey configurations, as labeled. Results are shown for the local (left panels) and orthogonal (right 
panels) bispectrum shapes (please note the differences in the vertical scales). The top two panels show the trend of the error as a function of & max (for fixed 
0.5 ^ z ^ 2.1 in the slitless case and 0.1 ^ z K 2.1 in the multi-slit case), the middle two as a function of the minimum redshift (for fixed & max = 0.3 Mpc~' 
and z ^ 2.1), and the bottom two as a function of the maximum redshift (for fixed k milx = 0.3 Mpc~' and z ^ 0.5 for the slitless case and z, ^ 0.1 for the 
multi-slit case). 



The improvement in the constraints on both /nl and cr 8 when 
combining the cluster and/or galaxy power spectra with the cluster- 
galaxy cross spectrum is reminescent of the fact that the latter is 
sensitive to different scales in a different way with respect to the 
former. The fact that this improvement is more important for the 
enfolded and orthogonal bispectrum shapes is extremely interest- 
ing, since these models are still relatively poorly studied compared 
to the equilateral and, almost ubiquitous, local shapes. Addition- 
ally, the amazing constraining power of the galaxy power spec- 



trum, even when considered alone, stresses the importance of the 
spectroscopy part for future Euclid-like missions when it comes to 
restrict cosmology on the basis of the spatial distribution of objects. 

The present work reinforces the conclusion that the spatial 
distribution of tracers of the LSS, especially galaxies, is an in- 
credibly powerful probe for primordial non-Gaussianity, thanks 
to the very strong impact t hat the latter has on the large scale 
bias of dark matter halos (jCarbone. Verde. & Matarrese 2008; 
ICarbone. Mena, & Verdell201oir The combination of cluster and/or 
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galaxy power spectra with the cross spectrum of clusters and galax- 
ies significantly improves the forecasted constraints, especially for 
the least studied non-Gaussian shapes. Merging all the information 
from future wide field surveys such as Euclid and WFIRST promise 
to bring constraints on /nl to the unity level, and constraints on cr 8 
to the l evel of ~ 1CT 4 , in both cases superior to future CMB exper- 
iments dSefusatti et al.ll2009l : IVerde & Matarresell2009h . 
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